Variational description of Mott insulators 
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The Gutzwiller wave function for a strongly correlated model can, if supplemented with a long- 
range Jastrow factor, provide a proper variational description of Mott insulators, so far unavailable. 
We demonstrate this concept in the prototypical one-dimensional t — t' Ifubbard model, where at 
half filling we reproduce all known phases, namely the ordinary Mott undimerized insulator with 
power- law spin correlations at small t' /t, the spin-gapped metal above a critical t' /t and small U , 
and the dimerized Mott insulator at large repulsion. 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.27.4-a, 71.30.+h 



Since Mott's original proposal the correlation- 
driven metal-Mott insulator transition (MIT) has at- 
tracted rising interest, renewed by the discovery of novel 
strongly correlated materials. On the verge of becom- 
ing Mott insulators, many systems display very unusual 
properties, high-T^ superconductivity being one spec- 
tacular example. While understanding Mott insulators 
and MITs is conceptually simple, calculations constitute 
a hard and long standing problem. Conventional elec- 
tronic structure methods, such as Hartree-Fock (HF) or 
density functional theory in the local density approxima- 
tion (LDA) cannot describe MITs, unless one allows for 
some kind of symmetry breaking. The standard example 
is long-range static magnetic order in the unrestricted 
HF, local-spin-density, or LDA-I-U approximations. This 
device works by effectively turning the MIT into a con- 
ventional metal-band insulator transition, thus masking 
the essence of the Mott phenomenon, where a charge gap 
appears quite independently of spin order. The fact that 
most known Mott insulators are indeed accompanied at 
low temperatures by some symmetry breaking, usually 
of magnetic type, further encourages the (wrong) sur- 
mise that it is not possible to describe any Mott insulator 
without a symmetry breaking. 

Another useful and popular approximation that may 
invite the same conclusion is based on the variational 
Gutzwiller wave function (GWF) and its various gener- 
alizations. H S S The GWF is the simplest way to 
improve a symmetry-unbroken, hence metallic. Slater de- 
terminant by partly projecting out the expensive double- 
occupancy charge configurations. In principle, were the 
projection complete, the GWF would indeed describe a 
Mott insulator devoid of symmetry breaking. Full pro- 
jection however means zero band-energy gain, generally 
incorrect, except for infinite on-site repulsion. For fi- 
nite projection, appropriate at finite repulsion, the GWF 
unfortunately describes a metallic state in any finite di- 
mensional lattice, at least so long as the uncorrelated 
Slater determinant state is metallic. 0] To obtain an in- 
sulator, one is forced once again to Gutzwiller project 



an artificially symmetry-broken determinant wave func- 
tion (WF). The main drawback of the GWF can be im- 
mediately recognized if one recalls Mott's original de- 
scription of a correlation-driven insulator. Let us con- 
sider for simplicity the single-band Hubbard model at 
half-filling. First of all it is clear that the system must 
even in the insulating phase allow for charge fluctuations 
around the mean value, implying some doubly-occupied 
(doublon, D) and empty (holon, H) sites, or it would 
jeopardize all band-energy gain. However, in order to 
describe an insulator, D and H have to be bound, other- 
wise any infinitesimal electric field could induce electric 
current. This spatial correlation among expensive charge 
configurations is exactly what is missing in the conven- 
tional GWF. There have been several attempts to cure 
this defects by allowing H-D correlations just on nearest 
neighbor sites, Q with only partial success. It is actu- 
ally suggested by our previous discussion that long-range 
H-D correlations should be a fundamental ingredient of 
a realistic insulating WF, otherwise there would always 
exist a finite probability for H and D to escape from one 
another, signaling of a metallic behavior. 

In this Letter we show that this plan can be accom- 
plished with success. In particular, we demonstrate for a 
specific one-dimensional (ID) model that: (i) non-trivial 
insulating WFs with or without symmetry breaking do 
exist; (n) the MIT, so far well described only in infinite 
dimensions, is accessible variationally even in finite 
dimensions. We consider the ID t — t' Hubbard model 

+ ^ X!"'*'T"»d. (1) 

i 

where (c^ ^) creates (destroys) an electron with spin 

a on site i and ni ^r — c|g.Cj ^. In the following we will 
assume t and t' positive. 

Bosonisation an d density-matrix renormalization 
group calculations predict that the ground state 

at half-filling is an insulator with gapless spin excitations 
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FIG. 1: Spin and charge properties of the COSO WF (top 
panels), the COSl WF (middle panels), and the CISO WF 
(bottom panels). Left panels: Nq/q (triangles) and Sq/q (cir- 
cles) as a function of q for L — 400. Right panels: D^, Zc, and 
Zs as a function of 1/L. The dotted line is a three parameter 
fit and the dashed line marks the value 1. The parameters 
of the WF are: Am = 0.33, A' = 0.17 {COSO and COSl), 
m = 0, A^ = 8.4 X 10"^ (CISO), Ai = 2t, Aa = 1.7t {CISO 
and COSO), and Ai = Aa = A3 = (COSl). 
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FIG. 2: Energy per site e for the simple GWF {L = 18, 
empty squares, and L = 82, full squares) and for the WF 
with long-range density-density Jastrow (GWF -I- J) (L = 18, 
empty circles, and L = 82, full circles), for t' = 0. The Slater 
determinant is a simple Fermi sea, with A^ — 0. The arrow 
indicates the exact energy per site of the Heisenberg model 
and the lines are guides to the eye. Top inset: accuracy of 
the WF with all the density-density Jastrow independently 
minimized (same symbols as before) and with the analytic 
parametrization (L = 18, empty triangles, and L = 82, full 
triangles). Bottom inset: Zc for the GWF -I- J (full circles) and 
for the GWF (open circles), lines are three parameter fits. 



(labeled COSl, where CnSm indicates a state with n 
gapless charge modes and m gapless spin modes) for 
t' /t ^ 0-5, a spin-gapped metal {CISO) with strong 
superconducting fluctuations for t'/t > 0.5 and small 
U/t, and a fully gapped spontaneously dimerized in- 
sulator {COSO) for f/t > 0.5 and large U/t. Dimer- 
ization is characterized by the discrete order parameter 
= lim|,_,|_^ Ix(«-i-l)-2x(i-j) + x(i-j + l)|, 
where x(i — j) ~ 9{Sf ■ Sf_^_i Sj ■ Sj_^_l) gives the dimer- 
dimer correlation function of spin rotationally invariant 
WFs, Sf being the spin operator along z-axis at site i. 

The CIS'O metallic phase suggests a variational WF 
built out of BCS, namely 



I*: 



= r\BCS)=Ve^p(J2f, 



|0), (2) 



where \BCS) is the ground state of a BCS Hamiltonian 
with gap function A, and dispersion = — 2cos(g) + 
2t' cos(2(7) — fJ,, fJ, being the free electron chemical poten- 



tial. The pairing function is defined by fq 
with the BCS energy spectrum Eg ~ 
A 



Ag/{eg+E, 



A2 

5' 



9^' 

and 



where Vn projects onto the subspace with fixed number 
of electrons N — L and J is the Jastrow factor: 



J = exp 



(3) 



where rii = Ui^^ + rii^i is the density operator, di = 
ni^lUi^l the doublon operator and hi = (1 — nij)(l — n^^i) 
the holon one. We remark that, in presence of a magnetic 
flux (j) ^ 0, the Jastrow term can acquire a non-centro- 



7^ 0. While the centro- 



- Wj i turns out to be less crucial 
ground state, the odd compo- 



symmetric component, 
symmetric part Wj j -|- 
than Vi_j for the < 
nent Wij — wj^i, which is zero at zero flux and had not 
been included so far, will play a relevant role in de- 
termining the response to (j). In practice we find it more 
convenient to discriminate the Mott insulator from the 
metal through correlation functions which can be calcu- 
lated in absence of fiux, as it is very unhandy to follow 
the changes of the variational WF due the presence of 4>. 
Numerically we find that the optimized Vij is smooth 



as a function of the distance 



while 



is stag- 



gered, Wij = \wij\ (— l)l'~-'~ I. Although Wij improves 
the accuracy, it does not change qualitatively the nature 
being variational parameters. Obviously \BCS) reduces of the WF, which is mainly determined by Vij, as we 
to the Fermi sea for A„ = 0. The operator P — JVn: discuss below. For this reason most of the results are 



g = Ai cos((j)-|-A2 cos(2q)-|-A3 cos(3g), Ai, A2, and A3 



3 



obtained with Wij = 0. Indeed, since (n; — l){nj — 1) = 
hihj + didj — hidj — dihj, then Vij < impHes H-H (D- 
D) long-range repulsion and H-D long-range attraction. 
The latter embodies the binding of H and D, while the 
repulsion prevents accumulation of H-D pairs. That is 
just the desired type of correlations missing in the GWF. 

We notice that optimization of the Jastrow factor al- 
lows to fulfill variationally the f-sum rule. Indeed, if we 
denote the variational energies of j^") by Eq and of Jigl^*) 
by £q , then we can prove for the optimized Vij 's that 



lim {Sq - £o) 
where iV„ = 



2iV„ 



(4) 



{^\n^qnq\^)/{'^\'^) is the equal-time 
density-density correlation function. On the assump- 
tion that our variational WF \'$) gives a fair descrip- 
tion of the actual ground state WF, Eq. Q allows to 
infer the behavior of the long-wavelength charge exci- 
tation spectrum by just analysing the charge structure 
factor Ng of j^*). Namely, if Nq ~ q, there should exist 
gapless charge modes, while if Nq ~ q^ charge excita- 
tions are presumably gapped. In turn Nq is determined 
by the structure factor of the uncorrelated BCS WF, 
N^, and by the Fourier transform of w^j , Vq, through 
Nq - N°/{1 + N°Vq). Q Since N^ tends to a constant 
value for g ^ 0, the long wavelength behavior of Nq 
is actually determined by the singular behavior of Vq, 
Vq ^ 1/q or Vq ~ implying Nq ^ q or Nq ^ q^, 

respectively. In other words the long-range Jastrow fac- 



tor not only restores gauge invariance, lim^^o Nq 



0, 



but also filters out either the metallic or the insulating 
component of the uncorrelated WF by an unbiased min- 
imization of the energy. 

Despite the large number of variational parameters, 
^i^j s, ^ijj s as well as A's, it is possible to optimize the 
full WF by recent developments in the variational Monte 
Carlo technique. |0| Yet, to afford larger scale simula- 
tions, we also assume an analytic parametrization of u^.j: 



A 



sin^(g/2) + — m 



(5) 



with A and m variational parameters, with m playing 
the role of the inverse Mott-localization length: m = 
or m > implying respectively metallic, Vq l/g, or 
insulating, Vq ~ l/g^, behavior. 

Besides the small-g behavior of the charge structure 
factor Nq, another quantity which can discriminate be- 
tween a metal and an insulator is a "rigidity" defined 
through Zc = (*|exp In the 

thermodynamic limit, for ID systems, \Zc\ — *■ for a 
metal and \Zc\ 1 for an insulator. Il6| 

For what concerns the spin modes, we argue that 
their behavior is determined by the BCS spectrum Eq, 
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FIG. 3: (a) \Zc\ as a function of U/t for t' /t = 0.75 and 
different sizes, (b) The same for the dimer order parameter 
Dd- The Vij's are independently minimized, very similar re- 
sults are obtained with the parametrization of Eq. @ . In the 
inset: the size scaling of \Zc\ for U/t = 5, the dashed line is a 
three parameter fit. 



namely if Eq is gapless then also the spin excitations 
are gapless and the spin structure factor Sq ~ q, 
whereas, if min(£'^) > 0, the spin excitation spectrum 
has a gap and Sq ^ q^. This is supported by the di- 
rect calculation of the generalized spin "rigidity" Zg ^ 

(*|exp(^E,j^|) |*)/(vl'|vl'), \Zs\ ^ for a spin- 

gapless WF and \Zs \ 1 for a spin-gapped WF. 

First of all we now prove that the WF (0) with ap- 
propriate choice of parameters is able to describe C05I, 
a Mott insulator with no symmetry breaking, as well as 
COS'O, a spontaneously dimerized insulator, and C'lSO, 
a spin-gapped metal. In Fig. ^ we show how the three 
possible cases COSl, (715*0, and (705*0 can be actually re- 
alized. The two WFs corresponding to (7050 and (7051 
are indeed insulating even though the average density of 
H (D) is finite. Moreover, we note that, unlike the un- 
correlated BCS WF, the correlated variational WF in the 
(7150 phase is not superconducting, as it should of course 
be in ID, thanks to the long-range Jastrow factor. \vi\ 
Remarkably, even if the (7050 BCS WF is not explicitly 
dimerized, after application of the Jastrow factor |(2Jl, a 
finite dimerization order parameter emerges. As in the 
liquid-solid transition of helium, it is typical of Jastrow 
functions to describe phase transitions in two-body corre- 
lation functions, without an explicit symmetry breaking 
in the one-body properties. 

To check the quality of the variational WF ^ for 
the Hamiltonian we start by considering the case 
of t' = 0, which allows us a direct comparison with ex- 
act Bethe ansatz results. For the reasons we mentioned 
previously, the GWF gives a rather poor variational de- 
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FIG. 4: Variational phase diagram of the ID t—t' Hubbard 
model. The error bars keep into account finite-size effects and 
dashed lines are guides to the eye. 



scription of the half-filled Hubbard model for U ^ t. In 
this limit an alternative procedure is commonly adopted. 
It is known that ior U ^ t the half-filled Hubbard model 
transforms by a Schrieffcr- Wolff transformation e' into 
an Heisenberg model. Therefore, instead of using the 
GWF, it is more convenient to use e*"^ \^g), where 
is the fully projected GWF. In that case, one searches 
for the optimal uncorrelated WF which, after complete 
projection, minimizes the average of e~*'^7ie''^, namely 
of an Heisenberg Hamiltonian. Practically one is forced 
to expand e*"^ ~ 1 + iS, which limits the applicability of 
this WF to the large U/t regime. 18] On the contrary our 
WF is equally accurate both at weak and at large U/t, as 
shown in Fig. |21 Remarkably the long range Jastrow cor- 
relations turn a fully metallic Fermi sea- the worst case 
one electron band for an insulator- into a perfectly good 
Mott insulator. 

Going back to the t ~ t' Hubbard model, in Fig. 13 
we show the behavior of \Zc\ and Dd for t'/t — 0.75. 
By increasing U/t, \Zc\ shows a sudden increase, and 
we can easily discriminate between the metallic phase, 
where \Zc\ ^ by increasing L, and an insulating phase, 
where \Zc\ — > 1 for L — > oo. Moreover, in the latter 
case, the MIT is indeed accompanied by a spontaneous 
dimerization, marked by a finite dimer order parameter 
Dd. see Fig.EIb). 

The variational phase diagram for the half-filled t—t' 
Hubbard model is finally shown in Fig. ^ In the region 
oit'/t < 0.5 and U > 0, we find no evidence of a phase 
transition, apart from finite-size effects at small U/t. The 
best variational state has A^, ~ g^, indicating a charge 
gap and Sq ~ q, gapless spin excitations: Aq connects 
only different sublattices (i.e., A2 = 0), making Eq gap- 
less. For t'/t > 0.5, there is a clear MIT at finite U/t 
between a spin-gapped metal, stable for small U/t and a 
dimerized insulator, stable at large U/t. In the metallic 
phase, the variational WF has Nq ~ q, Eq gapped, with 
A2 / 0, although |Ag| is small, hence Sq ~ g^, corre- 



sponding to an exponentially decaying spin-spin correla- 
tions. By increasing U /t and entering into the insulating 
phase there is a fast increase of A2, with Eq always fully 
gapped, and Vq ^ 1/q^. In this phase therefore Nq ^ q^ 
and Sq ^ q^. Although the variational WF does not 
explicitly break translational symmetry, there is a finite 
dimerization in the thermodynamic limit. Remarkably, 
this is due to the concomitant effect of a singular Jastrow 
Vq ^ 1/q^ and a gapped BCS spectrum Eq. 

In conclusion, it is possible to describe a Mott insula- 
tor by a variational WF based on a metallic Fermi sea 
which includes, besides a partial Gutzwiller projector, a 
long-range density-density Jastrow term. The accuracy 
of the WF is verified in a one-dimensional t — t' Hub- 
bard model at half-filling, which has a non trivial phase 
diagram comprising different phases and a MIT. We ar- 
gue that similar variational WFs with Vq ^ \/q^ could 
describe non-trivial Mott insulating phases in higher di- 
mensions too, giving the possibility to improve realis- 
tic descriptions of Mott insulating compounds. We also 
believe these wave functions to hold a great potential 
promise for future attacks to the superconducting phases 
that result by doping Mott insulators. 
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